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D ■ Abstract 

' ■ We discuss three complementary aspects of scalar curvature singularities: asymptotic 

' causal properties, asymptotic Ricci and Weyl curvature, and asymptotic spatial properties. 

We divide scalar curvature singularities into two classes: so-called asymptotically silent sin- 
gularities and singularities that break asymptotic silence. The emphasis in this paper is on 
the latter class which have not been previously discussed. We illustrate the above aspects 
and concepts by describing the singularities of a number of representative explicit perfect fluid 
solutions. 
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1 Introduction 



In a recent paper a general framework for locally analyzing Einstein's field equations (EFE) was 
presented [J . As a first application the structure of generic spacelike singularities was investigated, 
and an attractor describing the asymptotic dynamical properties of such singularities was presented. 
The results in P were later numerically supported in [2] and [3]. In ^ isotropic singularities 
were investigated by extending the methods developed in pp to also include the null geodesic 
\^ • equations. Gradually a picture about asymptotic causal structure has emerged, indicating that 

5^ ! generic spacelike singularities and isotropic singularities are examples in which local null cones 

collapse onto a timeline towards the singularity, leading to the formation of particle horizons 
with a size that shrinks to zero (see p. 75 and Fig. 1 on p. 76 in Since the shrinking of 

horizons prevents communication, the phenomenon was referred to as asymptotic silence and a 
corresponding singularity as being asymptotically silent (see |3] for further discussion). However, 
for some solutions particle horizons may not form, or if they do form, they may not shrink to 
zero size. The best known example showing that the particle horizon can be broken in one spatial 
direction is a spatially homogeneous (SH), Bianchi type I dust solution, with a so-called weak null 
singularity (see |^ p. 176). The asymptotic behaviour of this solution towards the singularity gives 
one illustration of how asymptotic silence-breaking can occur. 

Our goal in this paper is to illustrate the notion of asymptotic silence-breaking at a scalar 
curvature singularity by analyzing a selection of explicit perfect fluid solutions, thereby extending 
our understanding of the asymptotic causal structure of singularities. Although our examples 
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are non-generic, we believe that they may be of importance for more general phenomena, worthy 
of further study. In this context it is worth pointing out that non-generic solutions can play 
an important role in determining the behaviour of generic solutions. For example, self-similar 
cosmological solutions, which are certainly non-generic, have been shown to play an important role 
in determining the dynamics of generic solutions.^ 

The outline of the paper is as follows. In Sec. 2 we describe three complementary geometrical 
aspects of a matter singularity: firstly, according to its asymptotic causal properties in terms of 
asymptotic silence or asymptotic silence-breaking; secondly, according to its scalar Ricci and Weyl 
curvature; thirdly, according to how a spatial volume element is deformed towards the singularity. 
In Sec. 3 we give a number of selected examples that show how asymptotic silence can be broken 
and how this is connected with the above aspects. In Sec. 4 we conclude with some remarks about 
our results and their implications. A brief introduction to the orthonormal frame formalism and 
to Hubble-normalized variables is given in Appendix A. In Appendix B we give the non-tilted SH 
self-similar solutions, which describe the asymptotic behaviour of the examples in Sec. 3. 

Although we make use of the orthonormal frame formalism (see, for example jZj) for doing the 
necessary calculations, a detailed knowledge of this formalism is not required in the main body of 
the paper. 



2 Geometrical properties of singularities 

In this paper we consider perfect fluid solutions of EFE with a barotropic equation of state p = p{p), 
where p and p are the isotropic pressure and energy density in the rest frame of the fluid. Linear 
equations of state are characterized by the constant 7 defined by p = (7 — We assume that 
p > and < 7 < 2. 

The solutions that we consider have a matter singularity, i.e. the matter density p diverges 
since the overall length scale £ tends to zero 

lim p — +00 . (1) 

£~>0 

Our analysis of singularities relies heavily on the use of scale- invariant variables that are defined 
by normalizing with the appropriate power of the Hubble scalar H of the fluid congruence. The 
Hubble scalar is defined in terms of the 4- velocity u of the perfect fluid by 

H=lVau\ (2) 

and serves to deflne the length scale £ according to 

(3) 

where the overdot denotes differentiation along u. Important dimensionless quantities related to 
H are the deceleration parameter q and the spatial Hubble gradient r^, deflned by 

1=~Jp~^' ?-a = ^2 ■ (4) 

where ha'^ = da'^ + Ua projects into the 3-space orthogonal to u. 

We now introduce three complementary geometrical aspects of curvature singularities. 

^The Kasner solutions, which are self-similar, determine the past attractor for generic cosmological solutions 
(see 0), while the flat FL solution, which is also self-similar, underlies the generic phenomenon of intermediate 
isotropization (see 1111 . p. 312). 
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2.1 Asymptotic causal singularity structure 

The first aspect of singularities that we consider is the asymptotic causal properties towards the 
singularity. We define a part of a singularity to be asymptotically silent if all observers that 
approach it have particle horizons that shrink to zero size. If a particle horizon does not shrink 
to zero, or if a particle horizon does not form in one or more directions, we say that asymptotic 
silence-breaking occurs. 

Reference ^ heuristically links the formation of shrinking particle horizons to the limit 

hm Ea' = (5) 
(E) P- 103502-10) of the Hubble-normalized components of the spatial frame vector fields 

E.' = "4 (6) 

(these vector fields are introduced in Appendix A). However, as we will see, it seems to be necessary 
to also demand that 

hm f/a = , lim = , (7) 

in order to provide sufficient conditions for asymptotic silence to hold. Here the are the frame 
components of the spatial Hubble gradient and the Ua are the frame components of the 
Hubble-normalized acceleration of the fluid (see H117() and H121|) in the Appendix). 

Below we will show how asymptotic silence can be broken in various ways. We will give examples 
where 

lim. Ea = diag {const, 0,0), (8) 

where const ^ (if const appears in an asymptotic expression, it will from now on be assumed 
to be non-zero), leading to asymptotic silence-breaking in one direction only. We refer to such a 
singularity as partially silent. It is also possible that 

lim Eo,' = Ca' , (9) 

where the constant 3x3 matrix Cq.' has higher rank than 1 (rank 1 yields equation ijHJl), which leads 
to the breaking of asymptotic silence in more than one direction, although we shall not consider 
examples with this property, apart from the Friedmann-Lemaitre model to follow.^ In addition we 
give examples of asymptotic silence-breaking when lim^^o E^^ — oo for some values of a and i, 
and when some of the components of oscillate indefinitely as ^ ^ 0, so that lim^^o -Eq* does 
not exist. 

To illustrate some of the above let us consider the flat Friedmann-Lemaitre model: 

ds^ = -dt^ +t'^/^'' [dr^ +r^{de^ +sin'^9d(t>^)], (10) 
P = p=(7-l)p, u^d/dt. (11) 



The radial null geodesies are given by 



dt^ ^ 



^ =t4/37. (12) 
dr J 

^Asymptotic silence-breaking in more than one direction appears to necessitate violation of an energy condition, 
i.e. p < or even p -I- 3p < (see 1171 to follow). 
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Solving for r ~ r{t) shows that the past hght cone dX r — Q,t — \s generated by the foUowing 
family of null geodesies; 



, _ ^(37-2)/37 ^ _ .(37-2)/37 ,2 

r-b = I 37-2 ^ , " - 37-2 '7^3 (13) 

I — In t , 6 = In to ' 7 ~ 3 

Observe that 



\\m r = { , 3 ^ 2 ■ (14) 

t^o+ \ +00 , < 7 < f ^ ' 

It follows that if I < 7 < 2, the past light cone intersects the hypersurface i = at 

r = r,[t,)^^/^'^-'^'^\ (15) 

thereby creating a particle horizon, as shown in Fig. ^ Moreover, the horizon shrinks to zero size 
as the singularity is approached since limt^^o fh = 0, and thus asymptotic silence holds. However, 
if 7 < I no particle horizon is formed, see Fig. ^ and hence asymptotic silence is broken. 




Figure 1: The past light cone of the point P 
particle horizon at r = rh in the case 7 > |, 
7<|- 



and the initial singularity at t = 0, showing the 
and the absence of a particle horizon in the case 



In this example the Hubble-normalized spatial frame vectors are given by 

= diag (1,1,1), (16) 

so that 

r diag(0,0,0) , |<7<2 

^lim£;„^ = <^ diag (1,1,1) ,7 = 1 , (17) 

[ diag (+CX), +CX), +00) , 0<7<| 

which shows that in this example the condition (|SJ| is a reliable indicator of the existence of a 
particle horizon. This result is to be compared with the casual properties depicted in Fig. ^ 

2.2 Asymptotic scalar curvature singularity structure 

The second aspect of curvature singularities that is relevant to our discussion is the asymptotic 
behaviour of the Ricci and Weyl curvature scalars. 
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Firstly, it is of interest to consider the Hubble-normalized curvature scalars, the density pa- 
rameter 

^ = ^- (18) 
describing Hubble-normalized Ricci curvature (equivalently RabR""^ / H'^) , and 

r^ahcd r^^ r^ahcd 

i.e, the Hubble-normalized Weyl scalar and pseudoscalar invariants. If linif^^o 7^ 0, then the Ricci 
scalar curvature, or equivalently, the matter content, is dynamically significant at the singularity 
and if lim^^o C 7^ or lim^^o C* 7^ 0, then the Weyl scalar curvature is dynamically significant at 
the singularity. 

Secondly, it is of interest to consider the ratios 



/-^ r^ahcd /^abcd 

Cr .= -2 , .= -2 



(or equivalently CabcdC^"'^ / RabR"^ and C^^^^C'''"''^ / RabR""^) , indicating the relative rates at which 
the Ricci and Weyl curvature scalars diverge. Note that 

C C* 

= (21) 

as follows from p8|) - (|20(l . We here introduce the following classification. A scalar curvature 
singularity is said to be 

• Ricci-dominatcd if lim^^o Ct? = and lim^^o C/j = 

• Weyl-Ricci balanced if Cr , are bounded, and at least one of them has a non-zero limit, 
or oscillates indefinitely as £ — > 

• Weyl-dominated if one of Cr , C^ is unbounded, possibly oscillating 

(see IS] for a similar classification used for discussing future asymptotes, although C*ij^^C°''"^'^ was 
not considered in that paper). 

As an example, consider an isotropic singularity (see e.g., 0]). In this case the Hubble- 
normalized Ricci curvature is dynamically significant while the Hubble-normalized Weyl curvature 
is not, 

lim n = l, hm C = , lim C* = . (22) 

This leads to lim^^o Cr — and lim^^o C'^ = 0, and one thus has asymptotic Ricci dominance. 
On the other hand, at a generalized Mixmaster singularity (see ^), lim^^o = while the Weyl 
curvature is dynamically significant, which yields asymptotic Weyl dominance. Indeed, lim£_,oC 
and lim^_>oC* do not even exist, which implies that lim^^o C^r and lim^^g C^ do not exist either. 
The non-existence of these limits is due to the oscillating nature of the generalized Mixmaster 
dynamics. However, C and C* are bounded (since they are bounded on the Kasner subset and the 
vacuum Bianchi type II subset which describe the attractor for Mixmaster dynamics, see ITI) and 
hence C and C* exhibit bounded oscillations, while Cr and oscillate in an unbounded manner. 
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2.3 Asymptotic spatial structure 

In JI] (p. 30) a classification of matter singularities is given that reflects the change of shape of 
a spherical element of a fluid as a singularity is approached. This classification is based on the 
limiting behaviour of the scale factors ia in the eigendirections of the expansion tensor, and is 
applicable to models in which the limits of the scale factors exist or are infinite. In practice it is 
more convenient to base this classification on the Hubble- normalized expansion tensor (see e.g., 

m p- 19) 

00/3 = ^, (23) 
which is related to the Hubble-normalized shear tensor'^ Sq^ according to 

0a/3 = Sap + ■ (24) 

It follows that = 3. We are considering models for which the limit Oq,^ exists as the singularity 
is approached, say 

lhne„;3 = e^^, (25) 

where is a 3x3 symmetric matrix, which in general depends on spatial position. We diagonalize 
Q^fj, and let 8* denote its diagonal entries. If the limit of the deceleration parameter exists, 

hm q = q\ (26) 
then one can define the shape parameters pa according to 

3 3 



Po 



a=l '^-r 1 



which for the well-known Kasner solution become the so-called Kasner parameters. The exponents 
Pa are related to diagonalized limit values of T,af3 by 

The classification (first introduced by Thorne IF for Bianchi type I cosmologies) is based on 
the signs of the pa'. 

(i) (+,+,+), point, 

(ii) (+,+,0), barrel, 

(iii) (+,+,-), cigar, 

(iv) (+, 0, 0), pancake, 

and cycle on 1,2,3, cf. p. 121 ^J. Heuristically, if 6^ > (= 0, < 0, respectively), for example, 
the length scale in the 1-direction will tend to zero at the singularity (respectively a const, +oo), 
thereby justifying this terminology. 



2.4 Asymptotically self-similar singularities 

The simplest cosmological singularities are asymptotically self-similar, in the sense that the limits of 
all Hubble-normalized scalars exist as ^ — > 0, and equal the constant values of an exact self-similar 
SH solution, which we shall refer to as the asymptotic solution at the singularity. The examples 
in Section O except for the Szekeres and Wainwright-Marshman solutions in Sections 13.2.21 and 
13.2.31 have an asymptotically self-similar singularity. The corresponding asymptotic solutions are 
as follows: 



'^^ap is defined by 11171 in the Appendix. 
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i) the Taub form of flat spacetime ^2 P- 193], 

ii) the diagonal vacuum plane wave solution [ni P- 191], 

iii) the Collins VI/i perfect fluid solution p. 190]. 

One can think of the Taub form of flat spacetime as the special Kasner vacuum solution for which 
the curvature tensor is zero. For brevity we shall refer to this form of flat spacetime as the Taub 
solution. Likewise we shall refer to ii) as the plane wave solution. 

For any self-similar SH solution, the Hubblc-normalized expansion tensor Qap has constant 
components, and the deceleration parameter q is constant. These quantities determine the shape 
parameters Pa according to (|27|1 . For the three above-mentioned solutions we have: 

i) The Taub solution: 

= (1,0,0). (29) 

ii) The plane wave solution: 

(p„) = (1, r + y/r{l - r), r - y/r{l - r)), (30) 
where r satisfies < r < 1. The pa satisfy the conditions 

Pi+P2+P3=pI+pI+pI> I- (31) 

iii) The Collins VI^ solution: 

(p„) = (l,^^,to£), (32) 

where = (2 — 7)(37 — 2), with | < 7 < 2, and r satisfies < r < 1. The Pa satisfy the 
conditions 

^ =Pi+P2+P3>pI+pI+pI- (33) 

The line element for the Taub solution is given by 

^ ^ ^2^^2 ^ ^ ^^2 ^ ^34^ 

Full details about the solutions ii) and iii) are given in Appendix B. For a given cosmological 
solution with asymptotically self-similar singularity, one can identify the asymptotic solution by 
calculating the pa . 

3 Explicit examples 

We divide our examples into two main categories: spatially homogeneous examples and inhomoge- 
neous examples. The quantity C*i,c(fi°^^^'^ ^ ^^'^ ihvs C* and C]j, is identically zero for all examples 
except for the Wainwright-Marshman solutions in Section [3.2.31 Hence these quantities will not 
be given except in this last case. 

3.1 Spatially homogeneous examples 
3.1.1 Bianchi type I perfect fluid solution 

We consider the LRS Bianchi type I perfect fluid solution with a linear equation of state, given by 
([TT|, P- 199 with Pa = (1, 0, 0)): 

^ _^2(7-l) ^^2 ^ ^2 ^-l ^^2 ^ ^|(^y2 ^ ^^2) ^ (35) 
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where 

^2-7 =a + m^ t2-7 _ (36) 

The matter quantities are 

u = A-(^-i)a/at, p = P = (7-iV, (37) 

where a and m are two positive constants, and 7 satisfies < 7 < 2. 
The Hubble scalar is given by 

H = \A^-^t-^[l + A-'^^-'^mh^-^] ^t~^ ^00, (38) 

and the Weyl scalar by 

CafccrfC' ' = ^fjl^ ■ (39) 
There is a curvature singularity at t = 0, and the rates of growth of the curvature scalars are 

p ^ t-7 ^ 00 , CabcdC"'" - t-^'' +00 , (40) 

as follows from (|37|l and H^W|I . The shape parameters are Pa = (1,0,0), which by H29(l shows 
that the asymptotic solution at the singularity is the Taub solution. In addition it follows from 
Section 12.21 that the singularity is of the pancake type. Equations H38|l and (|40|l imply that the 
density parameter and the Hubble-normalized Weyl scalar satisfy 

n^t^-'^^O, C-t2(2-7)^o, (41) 

which is consistent with being asymptotic to the Taub solution. In addition, it follows from (|37|l 
and H39|l that the ratio Cr, as defined by (|20|l . satisfies 

(42) 

Equations H4()|l - (|42|l show that the singularity is Weyl-Ricci balanced, and is such that neither the 
Weyl nor the Ricci curvature is dynamically significant. 

It follows from H35|l and (|38|l that the Hubble-normalized frame components satisfy 

i;„'-diag(/3,t,i)^(/3,0,0), (43) 

3-1-2 

where P = 3a ^(^-t) ^ showing asymptotic silence-breaking in the x-direction, and this is further 
confirmed by the form of past null cone of a typical point, as shown in Fig. |21a). 

The dust case is discussed extensively by Hawking and Ellis on p. 145-147 in By means of 
a coordinate transformation that takes the Taub line element, —dt^ + dx^, to explicit flat form 
—dT^ + dX^, they show that one can make a C° extension of the metric across a null surface where 
the above "interior" solution is joined with flat spacetime (see our Fig. I^b and Fig. 22 on p. 146 
in |H]). For this reason the singularity is referred to as a "weak null singularity". 



3.1.2 Bianchi type VI/j tv^ro-fluid solution 

Wainwright ^5] has given an exact SH perfect fluid solution of Bianchi type Ylh with a barotropic 
equation of state. For our purposes it is useful to consider a special case of this solution determined 
by setting u = jq,v = w = ^ in |12| . in which case the group parameter is h — —1/9. In this 
case the matter content can be interpreted as two non-interacting dust and radiation fluids. The 
line-element is given by: 



ds^ = (-dt^ + dx^) + e^)'^/^ ^y2 _^ ^2 ea=)-2/5 ^^2 ^ 



(44) 
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Figure 2: Figures a) and b) show the causal structure associated with the LRS Bianchi I solution 
that exhibits a weak null singularity, and illustrates asymptotic silence-breaking. It is an example 
in which no particle horizon forms in one direction, i.e., it is an example of a partially silent 
singularity. Fig. utilizes the present coordinates while Fig. I^s) uses coordinates adapted 
to the null structure of the singularity. In both figures the thin dashed lines represent the fluid 
worldlines. 



where 

A = as+am t^l^ ■ (45) 
The constants a<,, (Xm are required to satisfy as > 0, OLm > 0. The matter quantities are given by 

u = ^ d/dt , PR + PM , P = ^PR , (46) 

where 

Pfl = i|«sa,„^-«/^ pM^^alt-^^^. (47) 

Note that asymptotically the equation of state is linear with a radiation equation of state, i.e. 
liint-^op/p =1/3. 

The Hubble scalar is given by 

H = lAr^ ^oo, (48) 

and the Weyl scalar by 

CabcdC^""' = -AH-^ [2(i)2r + OiT^)] , (49) 

where T — tdt\a A ^ t^/^ 0. There is a curvature singularity at t = 0, and the rates of growth 
of the curvature scalars are 

p ^ ^ oo , CabcdC^"'" - ~t-^^'^ ^ -oo , (50) 

as follows from H46|) . H47f) and (|49|l . The shape parameters are given by 

p„ = (l,|,-i), (51) 

and satisfy (|31(l . which shows that the asymptotic solution at the singularity is the plane wave 
solution with h = In addition, it follows from Section that the singularity is of the cigar 
type. 

Equations (|48|l and (|50|l imply that the density parameter and the Hubble-normalized Weyl 
scalar satisfy 
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which is consistent with being asymptotic to the plane wave solution. In addition, equation H5()|l 
shows that the ratio Cr satisfies 



-t 



-2/5 



-oo . 



which implies that the singularity is Weyl-dominated, in the terminology of Section 12.21 
It follows from H44|l and (|48|l that the Hubble-normalized frame components satisfy 

Eo,' ^ diag {const, t^'^) {const, 0, 0) , 



(53) 



(54) 



showing asymptotic silence-breaking in the x-direction. 

It is instructive to introduce the same "Taub" coordinate transformation as was used in [H] to 
show that the previous LRS Bianchi type I solution was C° extendible. Let 



T + X = < 
T -X = tc-- 



t = VT^ 



X2 



which leads to 



ds' 



A-^ {-dT^ + dX'^) + {T + Xf'^ di/ +A^{T + X)-2/5 



(55) 



(56) 



with A = as + am {T"^ — X'^Y/^ . The metric is C" on T - X = where it can be joined with flat 
spacetime in the region T ~ X < 0, but not on T + X = 0. Thus this plane wave singularity is 
also (partially) a weak null singularity. 




T+X=0 



T-X=0 



extension 



Figure 3: The causal structure associated with the Bianchi type VI^ two-fluid solution in T,X 
coordinates. 



3.2 Inhomogeneous examples 

3.2.1 Senovilla-Vera (1997) dust solution 

Senovilla and Vera have presented a number of explicit dust solutions in 13 . One which is useful 
for our purposes is case (i) on p. 3483 in |13| . This solution can be written in fluid comoving 
coordinates so that the line-element takes the form: 

ds^ = ^dt^ + <2 dx^ + y2 (fy2 _^ e-'^f-^ dz^ , (57) 

where 

Y = c^{t e-^'f- + {t e^f + c+ {t e-^ f+ , (58) 

b = i(l + /c), b± = i(l± V2~A:2), -1< A: < 1 . (59) 
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The parameters b, b± satisfy &_<0<&<l<fe+. The matter quantities are given by 

u = ,^ <'-^')";"-''' , P = 0. (60) 

The Hubble scalar is given by 

H = ^^[3-k + 2Y-^[c_b_{te-'')^- +m%{te^f + c+b+{te''^)''^]] , (61) 
and there is a simple expression for the Weyl scalar: 

2Y 



afccdL^ — P 



3 m? [t e^) 



(62) 



The solution yields three distinct types of singularity depending on the sign of c_. For our 
purposes, however, it suffices to consider the cases c_ > and c_ = 0. 



Case i) c_ > 

There is a curvature singularity at t = 0, and the rates of growth of the curvature scalars are 

P^t ^ — > OO , LabcdC —t ^ OO, [06) 

as follows from H60() and H62|l . The shape parameters are given by 

p„ = (l,fo-,l-&)- (64) 

It follows that the Pa satisfy H31|l . which shows that the asymptotic solution at the singularity is 
the diagonal plane wave solution, with h = —[1 — ^{k + y/2 — k^)]/[l + + y/2 — k'^)]. Since 
< < 6 < 1, the singularity is of the cigar type. 

It follows from Ht)l|l and H()3|) that the density parameter and the Hubble- normalized Weyl scalar 
both tend to zero: 

n^t''-''-^0, C - -t''-''- ^ , (65) 

as required for a solution that is asymptotic to the diagonal plane wave solution. In addition, 
equation H()3|) shows that the ratio Cr satisfies 

Cr ^ -t-^'-"-^ ^ ~oo , (66) 

which implies that the singularity is Weyl-dominated. 

This inhomogeneous solution is asymptotically homogeneous in the sense that the spatial Hub- 
ble gradient tends to zero: 

n - t''-''- . (67) 

It follows from H57|l and (|61|) that the Hubble-normalized frame components have the following 
limit: 

Ec' ~ diag A-^ , , t^^ ^ {const, 0, 0) , (68) 
which suggests asymptotic silence-breaking in the x-direction, confirmed in Fig. 0] 



3 EXPLICIT EXAMPLES 



12 



Case ii) c_ = 

There is a singularity at i = 0, with Y = also, and the rates of growth of the curvature scalars 

are 

p^t'^^oo, CabcdC'^''"' - -t-^ ^ -OO , (69) 
as follows from H6()(l and H62|l . The shape parameters are given by 

p„ = (l,i(l + A;),i(l-fc)). (70) 

It follows that the Pa are given by l|32|) with 7=1 and r = k, which shows that the asymptotic 
solution at the singularity is the Collins VI/i solution. Since — 1 < A; < 1, the singularity is of the 
anisotropic point type. 

It follows from Ht)l|l and H69|l that the density parameter and the Hubble-normalized Weyl scalar 
have the following limits: 

f]^|(i-fc2), c^-fii-e)\ (71) 

as required for a solution that is asymptotic to the Collins VI^ solution for dust (set j = 1, r — k 
in H143|) and H144|) ). In addition, equations (O and ij^U show that 

CR^-h (72) 

which implies that the singularity is Weyl-Ricci balanced. 
This solution is also asymptotically homogeneous, since 

ri^t''+-'' ^0. (73) 

It follows from H57|l and (|61() that the Hubble-normalized frame components have the following 
limit: 

E^' ^ diag (l , ^ ^(i+fc)/2^ ^ (^^^^^^ 0^ ^ (74^ 

which shows asymptotic silence-breaking in the x-direction. 

As in the SH cases it is instructive to make a "Taub" transformation of the coordinates 1551) . 
It turns out that the c_ > case is C° extendible through T + X = 0, but not through T — X — 
(the c- — case is not extendible) and thus the singularity in this case is (partially) a weak null 
singularity. See Fig. ^ for the causal structure of the singularity for the various cases. Asymptotic 
silence is broken in the x-direction and the singularity is partially silent for both cases. 

3.2.2 Szekeres solution 

The Szekeres solutions ^S] are a family of cosmological dust solutions. These solutions can 
be interpreted as exact perturbations of the Friedmann-Lemaitre models, with both a growing and 
decaying mode. For our purposes it suffices to consider a special Szekeres model whose line-element 
is 

ds^ = -dt^ + t^l'^ (X^dx^ + dy^ + dz^) , (75) 

where 

X^a + kxr'^, (76) 
and a, k are positive constants. The matter quantities are given by 



u = a/9t, p^\ar'^ x-^ , p = 



(77) 
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Figure 4: a) shows the causal structure associated with the SenoviUa-Vera 1997 solution in both 
cases c_ > and c_ = 0, in the original coordinates, while b) shows it for the c_ > case in T, X 
coordinates. 



The Hubble scalar is given by 



3t 



_^ kx 



(78) 



2tX_ 

and the Weyl scalar by 

CabcdC^"^' = ^p' ( 1 - - ) . (79) 



a 



If X < the singularity is determined hy X — 0, and the rates of growth of the curvature 
scalars are 

p^X-^^oo, C abcdC''''' X-^ ^oo, (80) 

as follows from H77|) and H79I) . The shape parameters are Pa = (1,0,0), which shows that the 
singularity is of the pancake type. 

It follows from H78|l and H80(l that the density parameter and the Hubble-normalized Weyl scalar 
satisfy 

rj-X^O, C-X^^O, (81) 
and in addition, equations 177|) and 179|l show that 

Cr^^. (82) 

which implies that the singularity is Weyl-Ricci balanced. 

The preceding results suggest that the solution is asymptotic to the Taub solution at the 
singularity. However, the behaviour of the spatial Hubble gradient shows that the solution is 
not asymptotically spatially homogeneous. It follows from H78|l and H121|l that there is one non-zero 
component given by 

= X{kx-2tXY ' 

which implies that 

n - x-i ^ oo . (84) 

It thus appears that the asymptotic state is not adequately described by the Taub solution. We 
shall refer to this singularity as being asymptotically "Taub-like" . 

It follows from H75|l and (|78|l that the Hubble-normalized frame components satisfy 

ii;„^^(/(x),0,0), (85) 
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where /(x) is a positive function, suggesting asymptotic silence-breaking in the x-direction. 

The causal structure in this example is quite different than in the previous examples. Even 
though Ei^ 7^ 0, a particle horizon does form for an observer with x < 0, and inititally shrinks 
in size as X ^ 0^ . However, the horizon does not shrink to zero size, as is confirmed by the 
numerical simulations that are shown in Fig. El Thus asymptotic silence is broken for the x < 
part of the singularity. 

This example also illustrates that the nature of the singularity can depend on the fluid worldline. 
If a; = a two-parameter set of fluid lines approach an asymptotically silent isotropic singularity 
(also known as an LK asymptote) [HI, 0] characterized by Pa = |{1, 1, 1) and O — > 1. If x > 
the fluid lines approach an asymptotically silent singularity of the cigar type described by the LRS 
Kasner asymptote = ^{—^, 2, 2) (see Appendix B). 




X 



Figure 5: This figure shows the particle horizon at an event P for an observer O, generated 
numerically. The particle horizon is determined by the past null geodesies through P, which 
intersect the singularity at a; = x^ and x = xj^. In case a) the particle horizon does not shrink to 
zero since as P approaches the singularity, x^ approaches Xcrit, the value determined by the dotted 
null geodesic. In case b) the particle horizon does shrink to zero since x^ xo as P approaches 
the singularity. 

As a result, the limits of Vi and C at the singularity are position-dependent (i.e. depend on 
x). In particular O ^ when x ^ Q but ^ 1 at a: = 0, i.e. as the singularity is approached 
develops a so-called spike (see also pp 91-96]). In addition, C ^ when a; < while C const 
when a: > 0, i.e., C asymptotically approaches a step function. 

3.2.3 Wainwright-Marshman solutions 

The stiff perfect fluid solution found by Wainwright and Marshman |18j can be written as 

ds^ = t^'^e''{-dt^ + dx^)+t^'^{dy + wdzf + f''^dz^ , (86) 
p = p= (m+^)t-2(™+i)e-n^ u = t-" e-"/2 , (87) 

Here w = w{t — x) \a an. arbitrary function, n = n[t — x) is determined according to 

n' = {w'f , (88) 

and m is a constant that satisfies m > — 
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The Hubble scalar is given by 

^^^-(m+l)g-n/2^^ (89) 

where 

S = i[(m+l) + it(u;')^]- (90) 
The non-zero Hubble-normalized components of the spatial frame vectors are 

The Hubble-normalized kinematic quantities are 

Finally the deceleration parameter and the non-zero component of the spatial Hubble gradient are 

The Weyl tensor can best be described by giving the complex Newman-Penrose scalars ^H] 

V-o = -imt-(™+i)e-", (94) 

^-4 = {~\m+\{w'ft + i{w" ~{w'f -mw' r^)t^l''\t-^'''^^^ e--^ (96) 
(see 122], p. 3037 for these formulae). The Weyl curvature scalars are then given by 

abcdC'^"^'^ - i C:,,<,C"''^'^ = 16 [i^l - 4V'i ^3 + V-o ^4) . (97) 
We choose the arbitrary function w{t — x) so that n(t ~ x) satisfies 

lim n{t — x) ~ — 00 . (98) 

The spacetime region is then defined by the inequalities 

t > , t-x>0. (99) 

Fluid worldlines with a; < encounter a curvature singularity at t = and those with a; > 0, at 
t ~ X, since p — + cx) in both cases, as follows from H87|l . In other words, the given initial singularity 
has two branches that depends on the sign of x, namely t = 0^ x < and t = x, x > 0, respectively. 
At the i = 0, a; < part of the singularity the solution satisfies 

lim ri = = lim Ui , lim (Ec-fl) = diag (2m - 1 , -i (4m + 1) , i (5 - 4m)) , (100) 

t->o+ t^o+ t~*o+ m -I- 1 ^ * 

which implies that the shape parameters Pa are given by 

Pa^^—im,\,l), (101) 
m -I- 1 * * 



and satisfy H131(l . The solution is thus asymptotic to a Jacobs stiff fluid solution. The singularity 

_3_ 

16 



type is: cigar (— -f^ < m < 0), barrel (m = 0), or anisotropic point (m > 0). We also have 



lim Ea' = 0, (102) 
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a) ^ b) ^ 

I I 



P 




singularity 



Figure 6: The figure shows the causal structure associated with the Wainwright-Marshman solu- 
tion. The left part of the singularity is asymptotically silent and is described asymptotically by 
the Jacobs stiff fluid solution. The right part of the singularity breaks asymptotic silence. 



and so we expect that the singularity is asymptotically silent, which is confirmed by Fig.|SJi). 

The nature of the t — x ^ 0,x > part of the singularity depends on the behaviour of the 
arbitrary function w(t — x) or n{t — x). However, since the null geodesies in the ix-space are given 
hy t±x — const, the causal structure in the te-space is independent of the nature of the behaviour 
of w(t — x) or n{t — x). Observe that the particle horizons of observers approaching the part of the 
singularity given by t — x = 0, x > do not shrink to zero size as the singularity is approached. 
Hence asymptotic silence is broken. Even though the causal structure is unaffected, the curvature 
properties depend significantly on the choice of w{t — x) and n{t — x). 

There are two possibilities that are compatible with equation H98|l : 

(i) limj^j,+ n'{t — x) = +oo 

(ii) n'{t — x) oscillates and is unbounded as i — > x"*". 
Case (i): \mit_,,j.+ n'{t — x) — +00 

It follows that 

lim B = +00 (recall n = (w'f) , (103) 

which in turn implies 

lim E„;3 = diag(2, -1,-1), lim g = 2 . (104) 

t—>x+ t^x+ 

It follows that the shape parameters are pa — (1,0,0), suggesting that the solution is asymptotic 
to the Taub solution. But equations H92|l and 1)93(1 imply that** 

lim ri = lim f/i = -3, (105) 

t—>x+ t^x+ 

*Here we make the additional assumption that Hmj_^^+ n"(t — x)/[n'{t — x)]^ = 0, which will be satisfied if, for 
example, n{t — x) = —C{t — x)~^, where b and C are positive constants. 
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so that the asymptotic state is inhomogeneous. It thus appears that the asymptotic state is not ad- 
equately described by the Taub solution. We shall refer to this singularity as being asymptotically 
"Taub-hke" (See Section E221. 

In addition, we have the following asymptotic expressions: 

H - e-"/2(u;')'-^oo, (106) 
p ^ e-"~>oo, CabcdC^'"' -^e-^^{w'f ^oo, (107) 

n - {w'y^^o, c-(w')"^^o, (108) 

Cr - {w'f ^ oo . (109) 

This part of the singularity is thus a Weyl-dominated scalar curvature singularity. Although it is 
"Taub-like" as regards the shape parameters, the curvature ratio Cji does not have the same limit 
at the singularity as in the solution that is asymptotic to the Taub solution (see equation (|42|l 'l. 
It is also worth noting that ^ 0, which shows that this solution is not in the generic class of 
stiff fluid solutions which satisfy 17 0, and whose behaviour near the singularity was analyzed 
by Andersson and Rendall 20 . 

Even though ^ as the singularity is approached, asymptotic silence is still broken, 
because the particle horizon to the left of an observer with a; > docs not shrink to zero at the 
singularity, i.e. does not tend to xq in Fig.|H|D). Thus this example shows that Ea^ —> does not 
necessarily lead to horizons that shrink to zero. This example thus forces us to make a distinction 
between Ea^ — > and particle horizons that shrink to zero. Since asymptotic silence should reflect 
an invariant physical causal property we therefore choose to define asymptotic silence as horizons 
shrinking to zero, cf. 3 . However, it should be pointed out that although the solution approaches 
the so-called silent boundary (see and 3,.), it approaches a part not described by SH equations 
(see and |21| for a discussion about the equations on the silent boundary) since f/a, 7^ 0. It 
is likely that this affects the relationship between E'^* and the causal structure. 



Case (ii): n'{t — x) oscillates and is unbounded as t x^ 
Let us consider the following possible choice |22| : 



w'it-x) = — ^ sin ( ), t-x>0 =^ (110) 

t — X \t — x J 

- -W^)^-b''-[T-x) ^'''^ 

(so that \mii^^+ n = —00, as required in equation (j^SJ). This yields 

• limt^2,+ Ea^ — {A, I ~ 2,3), while limt_^j.+ Ei^ does not exist, since Ei^ oscillate while 
remaining bounded. It follows that the singularity is non-silent. 

• Q, oscillates and is bounded. 

• limj__,^+ Yjap does not exist because of (bounded) oscillations. 

• C oscillates and is bounded but C* oscillates and is unbounded (due to w"). 

It is perhaps natural to refer to this kind of singularity as being strongly Weyl-dominated, i.e. 
the Hubble-normalized Weyl curvature is unbounded.^ This solution is thus a counterexample to 



^This behaviour at the singularity is reminiscent of the future evolution of Bianchi type VIIq models (see l'2'il ). 
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Solution 


asymptotic solution 


asymptotic spatial structure 


lim C R 


lim q 


-L-j-TlO iJldllClil Lypc 1 


TanVi 


pancake 


4 
3 


9 


Bianchi type Ylh 


plane wave {h = — i) 


cigar 


— OO 


3 
2 


two- fluid 










SV (1997) dust c_ > 


plane wave (— 1 < /i < 0) 


cigar 


—oo 


2 

l-3h 


SV (1997) dust c_ = 


Collins VI/i (-1 < /i< 0) 


point 


2 
3 


1 
2 


Szekeres 


Taub-like 


pancake 


4 
3 


2 



Table 1 : Summary of properties of the diagonal solutions that exhibit asymptotic silence-breaking 
at the singularity. 



the hope that the Hubble-normalized Weyl curvature is bounded at any cosmological singularity. 
It is nevertheless possible that the Hubble-normalized Weyl scalars are bounded at asymptotically 
silent singularities. 

4 Concluding remarks 

In this paper we have presented a collection of cosmological solutions that illustrate the phe- 
nomenon of asymptotic silence-breaking. The examples, apart from the Szekeres and Wainwright- 
Marshman solutions, are asymptotically self-similar at the singularity, in the sense that the limits 
of all Hubble-normalized variables exist as ^ — > 0, and equal the constant values of an exact 
self-similar SH solution. The properties of the singularities of these solutions are summarized in 
Tabled 

We comment on a number of common features of these examples: 

i) All of the solutions have a Ricci and a Weyl curvature scalar singularity in the sense that 

limp = oo, lim CabcdC'""^ ^±oo. (112) 

ii) Asymptotic silence-breaking occurs in one direction only. 
On the other hand, the differences are found in: 

i) The Hubble-normalized Weyl curvature components and scalars, and the density parameter: 

Taub £afi & Hafj ^0 C^O n^O 

Plane wave £af3 & Ti-ap t^O C— >0 fi— >0 

ColhnS Ylh £af3 & Ha/3 7^0 Cy^O VL^Q 

The variables £afj and T-Lafj are Hubble-normalized frame components of the electric and 
magnetic parts of the Weyl tensor as defined in Appendix A, equation (|118|l . 



ii) Asymptotic spatial structure (see Tabled) 



4 CONCLUDING REMARKS 



19 



iii) Weyl-Ricci balance, as described by the ratio Cu (see Table ^| 

It is of interest to ask what statements can be made in general about solutions that are asymp- 
totic to the Taub solution, the plane wave solution and the Collins VI/i solution at the singularity. 
On the basis of the properties of the asymptotic solutions in Appendix B, we can draw the following 
conclusions: 

• In any perfect fluid solution asymptotic to the Taub solution at the singularity, asymptotic 
silence-breaking occurs with rank (Ea^) — 1, for any 7, < 7 < 2. The spatial singularity 
type is pancake. We are unable to predict whether limf_,o Cfl = | general, since the ratio 
Cr is undefined for the Taub solution. A detailed analysis of the asymptotic behaviour would 
be needed to answer this question. 

• In any SH perfect fluid solution asymptotic to the plane wave solution at the singularity, if 
1 < 7 < 2 then asymptotic silence-breaking occurs with rank {E^) — 1, and the singularity 
is of the cigar type. The group parameter h will satisfy —1 < h < 0. We are unable to 
predict whether lim^^o Cr, — —00 ui general (if 7 > 1) since the ratio Cr is undeflned for 
the plane wave solution. A detailed analysis of the asymptotic behaviour would be needed 
to answer this question. 

• In any perfect fluid solution asymptotic to the Collins VI/i solution at the singularity, if 
1 < 7 < 2 then asymptotic silence-breaking occurs with rank {E^ ) — 1. The group param- 
eter h will satisfy — 1 < /i < 0. The spatial singularity type is point, barrel or cigar, and is 
necessarily a point if 7 = 1. We can also predict that if 7 = 1, the limit of the curvature 
ratio Cr will be — |. If 1 < 7 < 2, the limit will also be non-zero, but its value will depend 
on 7 and h through equations Ij 143(1 and H144|) . 

The remaining examples, the Szekeres solutions and the Wainwright-Marshman solutions (case 
ii)), are not asymptotically self-similar at the asymptotically silence-breaking part of the singularity. 
In the Szekeres solution, asymptotic self-similarity fails because lim^^Q+ ri = -l-oo, while the 
behaviour of the Hubble-normalized curvature scalars (51, C and Cr) is the same as the solutions 
that are asymptotic to the Taub solution. For this reason we refer to the Szekeres solution as being 
"Taub- like" at the singularity. In the Wainwright-Marshman solutions of class ii), the failure of 
asymptotic self-similarity is reflected in the curvature, in the sense that the limits of the Hubble- 
normalized energy density and Weyl tensor as ^ ^ do not exist. Indeed, C* is even unbounded. 
They also constitute an interesting example of asymptotic silence-breaking, since Ei^ oscillates so 
that the limit of this quantity does not exist. 

The Wainwright-Marshman solutions that belong to class (i) exhibit asymptotic silence-breaking 
even though E^ — *■ 0. In this solution, however, r^, Ua have non-zero limits. This raises the ques- 
tion: are the conditions ]iTai^o{Ea ,ra,Ua) = (0,0,0), given by equations ^ and Q, sufficient 
for asymptotic silence,^ i.e., the existence of particle horizons that shrink to zero as the singularity 
is approached? 
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®They are not necessary conditions since recurring spike formation (see pj) provides a counterexample; in this 
case, asymptotic silience holds, with lim^^g Ea^ = 0, but lim^^o and lim^^g do not exist. 
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A Orthonormal frame formalism and Hubble-normalized vari- 
ables 

In this paper we use the orthonormal frame formalism in conjunction with Hubble-normalized 
variables to calculate the properties of the explicit solutions. We introduce an orthonormal frame 
{eQ,ea}, where Eq = u, the 4-velocity of the perfect fluid, which is assumed to be irrotational. 
Local coordinates t, are introduced such that 

eo = N-^dt , e„ = ejd, , (113) 

where dt, di denote partial derivatives with respect to t and [a — 1,2,3, i = 1,2,3). 
The commutators are decomposed according to: 

[eo,e,](/) u^ea{f)-[H5j + aj + eJ-,W]ep{f) (114) 

[ea,e/3](/) = [2a^^5p{< + e^i,sn'-^)e^U) . (115) 

Here H is the Hubble scalar, u" the acceleration, aap the tracefree shear, 57" describes the angular 
velocity of the spatial frame { 6^ } relative to a frame that is Fermi-propagated along the integral 
curves of eo, while a" and n^^ = ?^(a/3) determine a spatial connection. Scale-invariant (dimen- 
sionless) Hubble-normalized variables are defined as follows (see PP; for a slightly different set of 
definitions based on conformal considerations, see |24| ) : 

do:=^=N-'dt, d^:=^ = E^^d,, Af:=NH, Ej:=^, (116) 

{[/", S^,^, A", N^p, i?"} {u-^, a^p, a", n^p, n"}/H . (117) 

It is convenient to express the Weyl tensor in terms of its electric and magnetic parts, relative to 
the timelike congruence Cq, denoted by Eap and Hap (see jllj . p. 19). Their Hubble-normalized 
components are defined by 

(EaPjHap) — {Eap, Hap)/{3H ) . (118) 

The Weyl scalars introduced in Section 12.21 are expressed in terms of Eap and Hap as follows: 

CabcdC^"'''' ^ HE^pE^^ ~ H^pH^^) , C:,^^C-'>-^ ^ l6E^pH"f . (119) 
It follows from H19() and H118() that the Hubble-normalized Weyl scalars are given by 

C = 72(£„^£"'3 _ -u^p-w-P) ^ c* = lU8^pn"f^ . (120) 

We note that the components Eap, Hap can be calculated using the orthonormal frame formulae 
in | 11| (p 35). Similar formulae for the Hubble-normalized components Sap, 'Hap are given in pp 
(p 21). 

In terms of the orthonormal frame formalism, the deceleration parameter q and the frame 
components of the spatial Hubble gradient Va, defined by are given by 

doH = -{q + l)H, daH = -raH, (121) 

It follows from Raychaudhuri's equation (see jlj) that 

g = 2S2 + i(37 -2)n-l {da -Ta + Ua- 2Aa) U°' , (122) 

where Y? = iSa^S"'^. The SH spacetimes are characterized by da{ ) = for all Hubble- 
normalized variables, and = = [/„, with the Cq, being tangential to the SH symmetry surfaces. 
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B The asymptotic solutions 

The asymptotic solutions that were introduced in Section lT^ are non-tilted SH self-similar solutions 
with perfect fluid source. A complete list of these solutions (line-element and matter variables) 
is given in jJT] (Section 9.1). In this Appendix we give the Hubble-normalized variables for the 
specific solutions that are used in this paper, which are the ones with diagonal line-element. Since 
the line-elements are diagonal, one can introduce an orthonormal frame in a natural way. For the 
line- element 

ds^ = -N^dt^ + X^dx^ + Y^dy^ + Z^dz^ (123) 

we choose 

eo = N-^dt, ei=X-'d^, e2=Y-^dy, = Z-^d, , (124) 
so that the non-zero components Cq* of the spatial frame vectors are 

ei^^X-\ e2^=Y-\ e3^^Z~\ (125) 
The self-similar solutions in question are either of Bianchi type I or type Vlh- 

B.l Bianchi type I solutions 

The hne-element is [m p. 187] 

ds^ = -df + t^P'dx^ + fP^dy^ + fP'dz^ . (126) 
The Hubble scalar and deceleration parameter are 

H^b-H-\ q = b-l, b=—^-^, (127) 

' ^ ' P1+P2+P3 ' ^ ' 

and the only non-zero Hubble-normalized commutation functions are the diagonal shear compo- 
nents, 

Sqq = bpa — 1 (no sum). (128) 
The constants Pa are the shape parameters (see ()28|) '). 

(1) Flat Friedmann-Lemaitre solution 

Pl=P2=P3 = ^- (129) 

The matter density is given by 

p=^t-\ 0<7<2, (130) 
which implies that the density parameter is = 1. The Weyl tensor Catcd is identically zero. 

(2) Kasner vacuum solutions and Jacobs stiff perfect fluid solutions 

^ ^ ^ 2_|^2^2l~lfoi' Kasner vacuum solutions, (131) 
Pi P2 P3 I Pi P2 ^■^ I fQp Jacobs stiff fluid solutions. 

The non-zero Hubble-normalized Weyl tensor components are 

£ii = '2p2 P3 - piip2 + Ps) ■ (132) 

where Sap — diag(£ii, £22, i^^as) and where £22,^33, are obtained from £11 by cyclically permuting 
1, 2, 3, and the density parameter is 

n = l{i-pl--pl-pl). (133) 

The choice (pa) = (1,0,0), and cycle, yields the Taub form of flat spacetime. 
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B.2 Diagonal Bianchi type VI solutions 

The line-element is p. 190] 

The Hubble scalar and deceleration parameter are 

H = b-H-\ q = b-l, 6=^^:^, (135) 
and the non-zero Hubble-normalized commutation functions are 

S„;3 = diag(6-l,6p2-l,6p3-l), ^1--|(C2 + C3), iV23 = -|(c2-C3). (136) 

The constants pi = 1, p2, P3 are thus the shape parameters (see (jSHJ)- The Hubble-normalized 
frame components are given by 

Ei^^b, ^bf-P^e""'''' , Es^ = bf~P^e-''''' . (137) 

The non-zero Hubble-normalized Weyl tensor components are 

£ai3 = dia.g{a,-^a + P,-^a- (3) , Uxi = ]^6^[(c2 + C3)(p2 -Ps) - (c2 - C3)(2 -Ps)] , (138) 

where 

a = ^b^ [(C2 - C^f - ip2 - Psf] , /3 = j^62 [id - 4) - (P2 - P3)(2 - P2 - P3)] ■ (139) 

(1) Collins Bianchi type Vlh perfect fluid solution 

2 — 7 ± rs (2 — 7)r ± s 



P2,3 



27 ' ' 27 

The parameter r determines the group parameter h according to 



C2,3 = ^ 7^ , s= v/(2-7)(37-2). (140) 



r^=rl{-h), re= \/(37-2)/(2~7), 0<r<l. (141) 
The matter quantities are given by 

!<7<S?. (142) 



which implies, using ()135|l . that 

o _ 

4 

The Hubble-normalized Weyl scalar is given by 



f^ = |(2-7)(l-r2). (143) 



C = -f 5^(1 - r') [(4 - 37)(2 - 7)(1 - r') + 127(7 " 1)] , (144) 

as follows from lfn5|) and ifT^ . 

The spatial singularity type depends on 7 and h. We restrict to non-negative pressure (7 > 1) 
which implies via p41|l that > 1 and hence that —h<r^. On writing the shape parameters p2 
and P3 in the form 

P2,3 = ^(l±rcr) (145) 
27 
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parameter range pa spatial singularity type restriction on 7 



rVc < 1 


(+, 


+,+) 


anisotropic point 


1 < 7 < 2 


rvc = 1 


(+, 


+ ,0) 


barrel 


1 < 7 < 2 


rVc > 1 


(+, 


+ ,-) 


cigar 


1 < 7 < 2 



Table 2: Classification of the singularity in the Collins Bianchi type VI/i perfect fluid solution with 
7 > 1. The value 7 = 1 corresponds to rj, = 1. 

using (|140|l and (|141|) . we arrive at the classification shown in Table |21 It follows from (|137|) that 
in all cases, 

lini Ea,' = diag(const, 0, 0). (146) 

(2) Diagonal vacuum plane wave solution 

P2,3 = C2,3 = r ± y/r{l-r) , (147) 
where r is related to the group parameter h through 

h = --^. (148) 
1 — r 

The Hubble-normalized quantities are given by 

Sii = y, , ^22,33 = feP2,3 - 1 , N23^-b^r{l-r), A, = -br. (149) 
1 + 2r 

The non-zero Weyl tensor components are 

£22 - -^33 = ^23 = 5&'(2r - l)^r{l - r) , (150) 

showing that the Weyl scalars p9|l . as given by p2()|l . are zero. 

The spatial singularity type depends on h (or equivalently r). The parameter r is in turn 
restricted by the value of the equation of state parameter 7 of any non-tilted perfect fluid SH 
solution that is past asymptotic to the plane wave solution, according to 

7<T^ (151) 

(see P- 663]). We restrict to 7 > 1, which implies that < r < ^ (since r > ^ implies 7 < 1). 
It follows from H147|) that P2 > and ps < 0, so that the singularity is of the cigar type. It also 
follows from H137() that 

lini Ec' = diag(const, 0, 0). (152) 
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